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Resonating valence bond states in a doped Mott insulator was proposed to explain superconduc- 
tivity in cuprates in January 1987 by Anderson. A challenging task then was proving existence of 
this unconventional mechanism and a wealth of possibilities, with a rigor acceptable in standard 
condensed matter physics, in a microscopic theory and develop suitable many body techniques. 
Shortly, a paper by Anderson, Zou and us (BZA) undertook this task and initiated a program. 
Three key papers that followed, shortly, essentially completed the program, as far as superconduc- 
tivity is concerned: i) a gauge theory approach by Anderson and us, that went beyond mean field 
theory ii) Kotliar's d-wave solution in BZA theory iii) improvement of a renormalized Hamiltonian 
in BZA theory, using a Gutzwiller approximation by Zhang, Gros, Rice and Shiba. In this article 
I shall focus on the merits of BZA and gauge theory papers. They turned out to be a foundation 
for subsequent developments dealing with more aspects that were unconventional - d-wave order 
parameter with nodal Bogoliubov quasi particles, Afileck-Marston's 7r-flux condensed spin liquid 
phase, unconventional spin-1 collective mode at (7r,7r), and other fascinating developments. Kivel- 
son, Rokhsar and Sethna's idea of holons and their bose condensation found expression in the slave 
boson formalism and lead to results similar to BZA program. 

At optimal doping, correlated electrons acquire a 'reasonable ' fermi sea, at the same time retain 
enough superexchange inherited from a Mott insulator parentage, ending tn a BCS like situation 
with superexchange as a glue ! Not surprisingly, mean field theory works well at optimal doping, 
even quantitatively. Further, t-J model is an excellent minimal model around optimal doping, where 
RVB superconductivity is also at its best. 
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Introduction 

Discovery of high Tc superconductivity in cuprates by 
Bednorz and MiiUer [1] in 1986 is a remarkable event. 
It was a breakthrough and a major turning point in the 
history of superconductivity and strongly correlated elec- 
tron systems. The field of quantum condensed matter 
physics and the community even got reorganized. About 
a month after Tanaka and collaborators confirmed[2] and 
brought to light Bednorz-Miiller's path breaking discov- 
ery, one of Anderson (PWA) jsj proposed a theory for 
high- Tc superconductivity. We followed it up with sev- 
eral papers [4, 5, 6,1 Si, focusing on superconductivity, 
within a span of one year. Several groups joined hands 
resulting in a flood of activities. 

Anderson's paper, to be referred to as PWA[3(], which 
suggested the resonating valence bond (RVB) mechanism 
of superconductivity, was very appealing and at the same 
time very unconventional. A challenge then was, as a 
first step, i) theoretically establishing RVB mechanism 
of superconductivity, in the t-J model and ii) develop a 
suitable many body theory that will be useful to calculate 
a variety of physical quantities and also confirm several 
of the predictions of PWA. This task was executed in 
two papers, in quick succession, by Anderson Zou and us 
(BZA)|3| and Anderson and us (BA)@|. During this pe- 
riod, except for a large value of superconducting Tc , not 
much was known experimentally, including fundamental 
properties such as, number of bands crossing the fermi 
level, symmetry of superconducting order parameter and 



survival of superexchange interaction in the doped Mott 
insulator. It is in this background, these two papers pro- 
vided a proof, acceptable in standard condensed matter 
physics, of the unconventional phenomena of RVB super- 
conductivity in a doped Mott insulator, and a theoretical 
frame work, a program, that became foundational for fur- 
ther quantitative developments. This was the beginning 
of the BZA program. 

The richness of Anderson's suggestion was that Super- 
conductivity in cuprates has turned out to be unconven- 
tional in more than one way: i) a new electronic mecha- 
nism, in an unexpected place called a doped Mott insu- 
lator and associated large Tc , ii) preformed spin pairing, 
later called 'spin gap' phenomenon, an unusual precur- 
sor to superconductivity, over a wide temperature region 
above Tc , iii) unconventional d-wave order parameter 
with nodal quasi particles, iv) unconventional spin-1 col- 
lective mode at (tt, tt) and v) an unusual competition 
from other types of charge and spin orderings etc. 

Superconductivity in cuprate family is a robust phe- 
nomenon at optimal doping. It has overcome disorder, 
charge and spin order tendencies and lattice instabilities. 
At optimal doping it is present in all cuprates containing 
Cu02 layers. A large condensation energy is evident in 
the way superconducting Tc , at optimal doping, jumped 
from the range of 30 to 90 to 120 and then 160 Kelvin, in 
new members of the cuprate family. A record Tc ^ 163^ 
is being held by a Tl based cuprate under a large ex- 
ternal pressure. As RVB theory was based on spins and 
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their exchange interactions, it was able to account for the 
large transition temperature and large condensation en- 
ergies in a natural fashion, compared to attempts based 
on phonon and other mechanisms. A strength of RVB 
theory from the beginning was its sound phenomenologi- 
cal basis, from where a flow of new concepts was natural. 
Mathematical difficulties that followed were in the na- 
ture of strongly correlated electrons in a tight binding 
band; a suitable many body theory did not exist. Inter- 
estingly, these formidable mathematical difficulties were 
also overcome, very efficiently, in the theoretical develop- 
ments that quickly followed. 

There have been efforts, before Bednorz-Miiller's dis- 
covery, in discussing possibility of superconductivity in 
models containing repulsive interactions, to understand 
superconductivity in heavy fermions and organics. His- 
torically, Hirsch[10| was the first to suggest an extended- 
s pairing in a repulsive Hubbard model. In a subse- 
quent paper Scalapino, Loh and Hirschllj interpreted 
the same superconductivity as spin fluctuation mediated 
pairing. Other authors have used the idea of d-wave pair- 
ing mediated by spin fluctuations in nearly antiferromag- 
netic metals, such as heavy fermions 12| and Bechgard 
salts [13|. The idea of pairing due to spin fluctuations 
continue to be pursued for high Tc cuprates by several 
One of the aims of the present paper is 
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groups 

to bring out, as clearly as possible, superexchange rather 
than exchange of spin fluctuations is a natural, physically 
correct and mathematically straight forward method to 
describe the glue for cuprates. Mott insulator is the 
template and superexchange is the glue. 

There has been a variety of efforts with varying suc- 
cess, in studying directly t-J and Hubbard models in 
2D, for superconductivity, along the RVB route: Kivel- 
son, Rokhsar and Sethna's ideaflGj] of soliton doping 
and bose condensation of holons, detailed slave boson 
analysis [7, 17, 3], detailed workfigj that sharpened the 
BZA phase diagram, an improved renormalized Hamilo- 
nian analysi s [201 . ga uge theory a ppr oaches [21. 
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32l . 13311 ■ variational monte- 
carlo l34ll . quantum monte carlo|35j, k-space[32ij and real 
space [37| cluster DMFT methods, diagrammatics 38, 
l39i powerful renormalization group studies _of_ Hubbard 

3, 



model in 2DWC 
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exact diaeonalization l42l |4S_ ^ 

series expansion [45|, |4^ and some analvtical|47l|methods 
have been employed. 

In the following sections, we elaborate a view that the 
basic and important problem of establishing an electronic 
mechanism of superconductivity by a many body the- 
ory and key physics of Mott insulators, during 1987, was 
solved in the very first phase. Anderson has expressed 
this view (isl . l49j in a recent article. The present article 
echoes similar views from a slightly different perspective 
giving some details. It has bit of history, as we will be 
completing 20 years since Bednorz-Miiller's path break- 
ing discovery. 



There is a recent review article by Lee, Na- 
gaosa and Wenjs^, which discusses physics of high 
Tc superconductivity in doped Mott insulators. It 
touches many of the theoretical developments. Our fo- 
cuss is on the BZA program, which has been so far very 
successful and has the potential to be used extensively 
for further quantitative progress. We also give some new 
insights and discussions. 

In the concluding section we contrast and distin- 
guish, theory of superconductivity in elemental metals 
from cuprates and other potential RVB superconductors. 
Even the type of questions raised and the way one ad- 
dresses issues are different. The scientific efforts put in 
unraveling the mystery of the complex cuprate system 
does not have many parallel in condensed matter physics. 
At the end, it is fair to say we do understand a lot, to 
be able to say where superconductivity comes from and 
what is the mechanism. There is more to be understood, 
of course. Such a realization has two effects: i) some what 
loose statement one hears occasionally, 'even now we do 
not understand high superconductivity and the mech- 
anism remains still unclear', looses its validity and ii) it 
gives one confidence and suggests that RVB theory is well 
and the BZA program is ready to answer new and old 
questions from experiments. Some encouraging recent 
examples from theory are: i) variational Monte Carlo 
analysis of Gutzwiller projected RVB wave functions by 
Nandini, Paramekanti and Randeria[34], and finding a 
good agreement with results of the BZA program as well 
as some experimental results ii) detailed calculation of 
electronic structure properties and excited state proper- 
ties by Gros-Muthukumar group|50], Ogata's group [51 
Zhang's group [Hll and others and iii) Anderson's very 
recent attempt to describe superconductivity and spin 
gap phenomena in an unified fashion using a notion of 
spin-charge locking and two types of Anderson-Nambu 
spinors. 

Brief Introduction to a Trio 



RVB proposal[3] in January 1987 created a spon- 
taneous involvement of theorists and experimentalists 
from all over the world, from widely different back- 
ground. The idea flourished instantly. Many key de- 
velopments took place during 1987-88. For example, 
Kivelson, Rokhsar and Sethna's idea[l3| of soliton doping 
and bose condensation of holons, spin-charge decoupling, 
Kotliar's d-wave solution [17|], Zou- Anderson's slave bo- 
son formalism^ adapted to t-J model, Affleck-Marston's 
phase 2l| , a condensate of 7r-flux of RVB gauge field, 
enlargement of U(l) to SU(2) RVB gauge theory[53|, 
Zhang- Rice singlet construction [HJl, an improvement of 
the renormalized Hamiltonian used in BZA theory [2^, 
sharpening of BZA phase diagram by detailed slave bo- 
son studi es (Til Il9l|. statistics transmutation[55l| . Laugh- 
lin's idea[56| of semion superconductivity (condensation 



of holons carrying an RVB gauge field flux of -j), other 
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Fukuyama, .. 

Slave Boson 
mean field theory 



FIG. 1: Evolution of BZA program beginning with Ander- 
son's proposal. Kivelson et al.'s bose condensation of holon 
is an independent development that also contributed to BZA 
program through slave boson approach. 
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TABLE L Comparison between conventional BCS theory and 
1987 RVB theory. BCS = Bardeen-Cooper-Schrieffer theory, 
RVB — Resonating Valence Bond theory, QP — quasi parti- 
cles, GL = Ginzburg-Landau theory, BdG = Bogoliubov de 
Gennes Equations, MFT = mean-field-theory. 



These intense activities revolved around PWA, BZA 
and BA papers. We call these 'trio', as these three pa- 
pers have a close link and continuity. The RVB super- 
conducting state and related theoretical developments 
in the trio in 1987 is similar to a BCS type of the- 
ory, but based on electron correlation mechanism, with 
its own novel features and notions and some formidable 
theoretical problems. This is shown in Table 1. After 
phonon pairing mechanism and BCS theory, RVB the- 
ory is the most significant development in the field of 
superconductivity, involving an entirely different mech- 
anism based on electron correlations and more impor- 
tantly compelled by a rich phenomenology. Even though 
heavy fermion superconductors and organic superconduc- 
tors existed around 1987, with an electronic mechanism 
at work, it did not excite the condensed matter commu- 
nity as much as cuprates did. 

1. Anderson's Proposal 

A qualitative and quantitative understanding of high 
Tc superconductivity in cuprates involved, first identi- 
fying the predominant mechanism of superconductivity. 
This, in turn, involved three major steps: i) abstracting 
key notions and introducing a new paradigm, ii) iden- 
tifying the right low energy effective Hamiltonian and 
iii) developing suitable theoretical methods and finding 
approximate solutions. Finding a new paradigm, ab- 
stracting new notions and a model for cuprate supercon- 
ductivity, using known phenomenology of La2Cu04 and 
other magnetic oxides is a remarkable chapter in con- 
densed matter physics. As PWA has expressed in an 
articlejll] entitled 'Magnetician's edge', a detailed and 
in depth knowledge of quantum magnetism (particularly 
in oxides) was essential. The phenomena of cuprate su- 
perconductivity turned out to be a meeting ground of 
quantum magnetism and superconductivity. The spirit 
of the approach to this complex quantum condensed mat- 
ter problem is well summarized in 'Central dogma in high 
Tc superconductivity', a chapter in PWA's book on high 
Tc superconductivity [9|. 

PWA^ identified the parent compound La2Cu04 as 
a spin-i Mott insulator, having one electron in a non- 
degenerate orbital per copper atom. This as well as the 
doped La2_xBaxCu02 is described by a single band large 
U Hubbard model in 2D at and away from half filling: 



h.c. 



(1) 



PT violating chiral RVB states [57[ Chern-Simons gauge 
fields, computation of physical quantities using RVB 
gauge field theory, anomalous normal state, failure of 
fermi liquid theory, electron confinement, inter layer pair 
tunneling and more. Some of them that are directly re- 
lated to one layer superconductivity are shown in figure 
1. 



Here the site index refers to a Wannier orbital. It is a 
symmetry adapted hybrid of Cu 3(ia:2_y2 and oxygen 2p 
orbitals, that retains d^^-yi symmetry. At half filling 
and when U >> t, the ground state is a Mott insulator. 
It has a finite Mott Hubbard gap for charge carrying 
excitations. Various estimates of U and t and also next 
nearest neighbor hopping's exist now: U ~ 5 eV, t ^ 0.25 
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eV. The ground state of the Hubbar model at half filling 
is well approximated by the ground state for hopping 
matrix element t = : 

|ai,a2,...,a^)^ J] c^jO) (2) 

i=ltoN 

In these states, every site is singly occupied and has a 
dangling spin. Consequently, total spin degeneracy of 
this manifold is 2^. The extensive spin entropy of the 
above states are removed by superexchange, a second or- 
der hopping processes, involving two neighboring sites at 
a time. By a second order perturbation procedure we 
can derive an effective Hamiltonian that lifts the 2''^ fold 
spin degeneracy. For a given pair of neighboring sites, 
the ground states for hopping, t = 0, are: 

IT,i), U,T), IT,T), U,i), (3) 

These 2^ neutral spin configurations are three bond 
triplets and a bond singlet state of two spins of neigh- 
boring sites. When hopping t is introduced perturba- 
tively, there is a virtual transition or mixing of the above 
states with the excited 'ionic spin singlet' intermediate 
configurations: 

I n,0} and |0,n), (4) 

resulting in a bond singlet ground state and bond triplet 
excited state. This is represented by the following effec- 
tive Hamiltonian for the spin dynamics of the large U 
repulsive Hubbard model: 

ir(half filling) ^Hs = J^(S, • S, - i) (5) 

(y> 

where J = 

PWA suggested that superexchange survives in the 
doped Mott insulators, up to some reasonable dopings. 
Electrons delocalize, but continue to respect the double 
occupancy constraint. The resulting model for non-half 
filled case is t-J model that contains in addition to the 
superexchange term Hs also the kinetic energy term H^: 

Htj = Ht + Hs (6) 

{ij) («i> 

with a double occupancy constraint, n^j + n^j^ ^ 2 at 
every site. Cuprates have a large superexchange J ~ 
0.15 ev, one of the largest among spin-i Mott insulators. 

PWA further suggested, in view of strong quantum 
fluctuations arising from spin-^ character and 2 dimen- 
sionality, that ground state of this 2D Heisenberg model 
is a quantum spin liquid, with a possible psuedo fermi 



surface for certain neutral fermion excitations. The mag- 
netic susceptibility data of Ganguly and RaojH^], for in- 
sulating La2Cu04 , which did not exhibit any phase an- 
tiferromagnetic phase transition feature, also seemed to 
support PWA's earlier notion of spin liquids 60] in 2D 
spin-i Heisenberg antiferromagnets. 

Singlet correlations in this quantum spin liquid was 
suggested as the neutral singlets or preformed pairs that 
are waiting to superconduct, given an opportunity. On 
doping, a fraction x of neutral resonating singlets get 
charged resulting in superconductivity. Here x is the 
doping fraction. The RVB mechanism was expressed suc- 
cinctly in the form of a Gutzwiller projected (double oc- 
cupancy removed) BCS type wave function 

|RVB;(/)> = PGUiuk+Vkcld.^m 

k 

that nicely interpolates the spin liquid ground state of the 
Mott insulator (x = 0) and the superconducting doped 
Mott insulator (x 7^ 0). Here N is the number of sites and 
N(l-x) is the number of electrons. Gutzwiller projection 
is defined as, 

Pg - n,^n,^) (8) 

i 

There was another important suggestion in this paper: 
two electrons in a given cooper pair will avoid double oc- 
cupancy and cooper pair function = for i = j. 
This automatically allowed extended-s and higher angu- 
lar momentum symmetry such as d-wave. 

This paper was very special. Superconductivity 
emerged from a Mott insulator (non-fermi liquid) back- 
ground. Pairing was not in k-space: superexchange, an 
intrinsically real space quantum chemical binding phe- 
nomenon lead to zero momentum condensation of cooper 
pairs. This paperfs'l has become a classic in quantum 
condensed matter physics, almost a poem that opened 
a new door and one that reveals new shades of meaning 
each time one reads it. 

2. BZA theory 

BZA followed heels and provided a physically motivated 
approximation method for the strongly correlated elec- 
tron problem at hand. It ended up being a begin- 
ning of a program. As mentioned earlier, PWA sug- 
gested a Gutzwiller projected variational wave functions 
parametrized by a pair function (t>(ij)- This theory un- 
dertook this variational analysis. This is similar to 
a BCS-Hartree Fock type analysis, but in a restricted 
Hilbert space containing no double occupancy. That is, 
one would like to minimize the energy expectation value 
(or free energy) with respect to the pair function (t>{ij): 

E[4>] = (RVB; ^{PaiHt + H,)Pg\RYB; 4>) (9) 
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Presence of Gutzwiller projector Pq rnakes computation 
formidable. This theory introduced a physically moti- 
vated approximation. The approximation amounts to 
treating the Gutzwiller projection in a mean field fashion 
and approximate the above expression by 



E[(P] « {RVB; (P\{xHt + Hs)\RVB; 



(10) 



That is, the complicated Gutzwiller projection was ap- 
proximated by replacing the hopping parameter t by a 
renormalized parameter xt, since x is the probability that 
an electron can find a neighboring site empty to which it 
can hop. In other words we have a renormalized Hamil- 
tonian 



Htj = xHt 



(11) 



defined in the full Hilbert space, also containing dou- 
ble occupancies. The rest is very similar to standard 
BCS theory. Interestingly this paper conjectured that 
this renormalization prescription should work well be- 
yond about 5% doping, about which we will discuss later. 

Within the above mentioned approximation this theory 
found i) a spin liquid ground state, neutral fermion exci- 
tations with a pseudo fermi surface for the Mott insulator 
and ii) a superconducting ground state with extended-s 
symmetry for doped Mott insulator. 

A key point in BZA paper is a liberation from the Pauli 
spin operators, that is traditionally used for analysis of 
quantum spin systems and go to the constituent electron 
variables, even for the Mott insulator, by rewriting the 
Heisenberg part of the Hamiltonian as 



(12) 



with the pair potential having the form, 7(k — k') ~ 
[cos{kx — k'x) + cos(fcy — k'y)]. It should be noted that the 
pair potential, while it is attractive for small momentum 
transfer (k — k') ^0 , changes sign and becomes repul- 
sive for large momentum transfer ~ (7r,7r), manifestly 
suggesting a dx2_j,2-wave rather than extended-s wave 
as a low energy mean field solution. For a reason that 
will be elaborated later, we were extremely satisfied with 
extended-s wave mean field solution. 

It then employed a Bogoliubov-Hartree-Fock factoriza- 
tion and identified nearest neighbor self consistent pa- 
rameters: 

A = ^(cosfta; -I- cosfcy)(c^|clj,^) and 

k 

P = "^((^oakx + coaky){cl^Ck„) (15) 

k(T 

The first one A is the usual anomalous superconducting 
amplitude. The second one p is somewhat unconven- 
tional, it is a kinetic energy or hopping term, a 'Hartree- 
Fock Vector Potential', generated by superexchange pro- 
cess. This unusual Hartree-Fock factorization term in- 
troduced in this paper played crucial role in later devel- 
opments, such as gauge theory^5] and AfFleck-Marston's 
flux phase [21 1 . 



2.1 Mott insulator 

Let us first consider the case of zero doping, x = 0, the 
Mott insulator. The simplest self consistent solution was 
found to be A = 1 and p = 0. This resulted in a quasi 
particle Hamiltonian for neutral fermions, a's : 



(ij) 



niF 



J I cos 



COS k,, 



(16) 



using the relation, Si = J2a,p4a'^a(3C^i3, where t is the 
Pauh spin operator. Further = ■^(Cit'^ii ~ ''4'^]t) 
the bond singlet or (in the present case) neutral cooper 
pair operator. In the electron representation the Heisen- 
berg Hamiltonian has a simple meaning. The spin-spin 
coupling encourages bond singlets, because it is minus of 
the bond singlet number operator bl^b^y The non-trivial 
character of the lattice problem arises from the fact that 
the bond singlet number operators do not commute, if 
they share one common site. We [gH showed recently a 
very useful commutation relation. 



(13) 



. It has a deep meaning that singlet resonance or delocal- 
ization involves an unavoidable spin chirality fluctuation. 

In k-space the cooper pair scattering term arising from 
superexchange has the following form: 



(14) 



The pseudo fermi surface for the neutral fermions is given 
by the expression \coskx -t-cosfcj,| — 0. Further, the 
anomalous pairing leads to a remarkable result for ground 
state occupancy 



(17) 



Even though neutral fermion excitations have a pseudo 
fermi surface, there is no momentum space discontinuity 
for the constituent electrons. In this sense this spin liquid 
ground state of the Mott insulator is far removed from 
any standard fermi liquid state. The ground state sug- 
gested by this RVB mean field theory is the Gutzwiller 
projected spin liquid state: 



I2DRVB) =PGn( 



Uk 



Vkcl^cl 



)|0) (18) 



k,k' 



with = ±1 inside and outside the pseudo fermi sur- 
face respectively. Emergence of neutral fermions with a 
pseudo fermi surface, in a many body theory for quan- 
tum spins, was a great excitement at that time. It was 
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radically different from bosonic spin wave excitations in 
ordered antiferromagnets. 

This paper also pointed out that in the Mott in- 
sulator, 'neutral fermion' quasi particles are meaning- 
ful only when they are created as 'particle-hole' pairs 

nited on opposite side of the fermi surface. When they 
are on the same side of the fermi surface we end up cre- 
ating charged fermions, which are not part of the low 
energy excitation spectrum of the Mott insulator. This 
implied that Brilluoiun Zone of neutral fermions are only 
half of the full BZ, very much like the BZ of spinous in 
the case of ID Heisenberg spin-i antiferromagnet. 

2.2 RVB superconductor 

The mean field analysis was then performed for the 
doped Mott insulator, x 7^ 0. In addition to the superex- 
change term it had the renormalized kinetic energy term 
xHt (equation 11). The problem at this level becomes 
very much like a standard BCS analysis, with a renor- 
malized band width and nearest neighbor cooper pairing 
of strength J. Thus superexchange becomes the glue. 

In addition to the anomalous average A 7^ (unlike 
the case of Mott insulator) it also found p / 0. It is 
a renormalization of kinetic energy from superexchange 
term. The overall superconducting solution corresponds 
to a spin singlet superconducting state with extended-S 
symmetry. 

The mean field solution suggested superconductiv- 
ity (ODLRO) for any filling. However, at half filling, 
ODLRO is destroyed by a complete suppression of low 
energy charge fluctuations, implemented by Gutzwiller 
projection. That is, at half filling the mean field super- 
conducting state turns into a neutral spin liquid state 
upon elimination of charge fluctuations. The mean field 
solution suggested superconductivity (ODLRO) for any 
filling. However, at half filling, ODLRO is destroyed by a 
complete suppression of low energy charge fluctuations, 
implemented by Gutzwiller projection. That is, at half 
filling the mean field superconducting state turns into a 
neutral spin liquid state upon elimination of charge fluc- 
tuations. 

It was pointed out that the mean field transition tem- 
perature obtained by this analysis, which has the generic 
BCS form, 

ksT: - J (19) 

a fraction of the superexchange J, is the cross over tem- 
perature at which 'preformed spin singlet pairs' are be- 
ginning to appear in the system. In current terminology 
this temperature scale is the 'spin gap' scale, below which 
spins are progressively getting paired in a cooperative 
fashion. These singlet pairs start transporting charge 
and compete with single electron transport. Thus the 
mean field temperature provides some kind of umbrella 



below which the charged valence bonds can undergo BEG 
type of condensation on their own, making use of their 
light mass and small density. For small doping the actual 
superconducting transition temperature will be bounded 
by the mean field Tc . 

As doping x increases, the effect of superexchange is de- 
creasing and at the same time a fermi sea is being built, 
because of increase in electron delocalization. Thus we 
have some kind of fermi sea with superexchange as the 
glue. So the Tc predicted by mean field theory will be- 
come close to the superconducting Tc . This is shown 
in figure 2, which reproduces the BZA mean field phase 
diagram. Beyond the dashed line, mean field results for 
superconducting Tc starts making sense. It is interest- 
ing to see that this is close to optimal doping, that was 
discovered in experiments months later. 
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FIG. 2: Phase diagram reproduced from BZA 4] paper. This 
phase diagram, drawn in March 1987, was based on a micro- 
scopic calculation plus some conjecture about validity of mean 
field theory at higher dopings. Surprisingly, the conjectured 
region of validity of meanfield theory (beyond dashed line) 
very nearly corresponds to optimal doping of ~ 15 % and 
beyond, that was experimentally confirmed later. 

At small doping superconductivity was viewed as a 
condensation of a fraction x of valence bond bosons that 
are charged and gave a BEG type of expression of Tc for 
small doping: 

Here m* is the effective mass of the charged valence bond 
and vq is the volume occupied by each Cu cell. What is 
remarkable about this formula is that the superconduct- 
ing Tc has a strong and explicit dependence only on t and 
x. It does not have any explicit dependence on J or the 
Hubbard U !. The large U has done its job through su- 
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perexchange, of preparing spin singlets at a sufficiently 
high temperatures. When BEC takes place below spin 
gap scale keT*, spin singlet correlation is maximal, and 
it is a coherent charge fluid, without much spin activity 
at low energies. This formula was modified, to a more 
appropriate Kosterlitz-Thouless type formula in 2D, in a 
subsequent paper Q by with Anderson, Hsu and Zou, as 

ksT, -ix-xc) 21 

m* 

Here Xc is some critical doping needed to overcome disor- 
der effects and begin superconductivity; this expression 
was taken from a similar expression for superfluidity of 
''He in vicor glass. At very high doping Mott insulator 
turns into a (disturbed) fermi sea. Superexchange be- 
comes less relevant, as electrons are less localized because 
of decreasing correlations. This mean field theory showed 
a sharp decrease of the mean field beyond an optimal 
doping. Synthesizing various ideas and the BZA mean 
field solutions, aphase diagram was suggested shortly, in 
the same paper [6|. This phase diagram, shown in figure 
3, was also a prediction of BZA theory. The experimen- 
tal phase diagram that was established later over years 
has a striking resemblance to this prediction. 



RVB meon field transition 




x=002 x=020 

Doping fraction (X) 



FIG. 3: Phase diagram reproduced from a paper by Ander- 
son, GB, Hsu and ZoujGj). This phase diagram, drawn in April 
1987, was based on a microscopic calculation plus some con- 
jecture, before a systematic study of doping versus Tc began. 
The overall agreement with the current cuprate phase dia- 
gram is striking. 

This idea of bose condensation of charge valence bonds 
was given a sharper expression with important conse- 
quences, as a holon (soliton) condensation by Kivelson, 
Rokhsar and Sethna[16|] shortly. 

The importance of quantum phase fiuctuations in a 
condensate of preformed pairs was emphasized in the 



BZA paper; these ideas/predictions which predated de- 
tailed experiments, have matured into notions such as 
phase fiuctuation dominated bad metal phase [il] and 
vortex liquid phase f63l|. with help from a variety of ex- 
periments performed after 1987. This basic ideas from 
RVB theory was used in an important paper by Emery 
and Kivelson (g^ to talk about quantitative comparision 
of superconducting Tc with experimentally measurable 
quantities, such as superfluid stiffness, using a phase fluc- 
tuation dominated scenerio. We will discuss this later. 

It is also worth emphasizing that BZA paper consti- 
tutes an unequivocal prediction of the spin pseudo gap 
phenomenon, long before it was noticed experimentally. 

Even though mean field theory was in agreement with 
experiments in giving high superconductivity at op- 
timal doping, an important question was whether the 
mean field superconductivity will survive quantum fluc- 
tuations arising from double occupancy constraint, this 
paper conjectured that when the doping is above w 5% 
the fluctuation correction arising from the double occu- 
pancy constrant should not matter very much. However 
it was necessary to validate this in an acceptable fashion 
either using existing many body theory or some thing 
new. 

3. Gauge Theory 

After BZA theory one of the urgent job was to give 
mathematical expression to the increasing phase fluctu- 
ations of cooper pairs, as one approaches the Mott insu- 
lator by decreasing the doping. It also amounts to find- 
ing how to express mathematically difference between a 
Mott insulator and doped Mott insulator. The gauge 
theory paper Q undertook this study from a very differ- 
ent perspective. It offered a new approach, which dis- 
tinguished Mott insulator and doped Mott insulator and 
enabled a systematic study of Mott insulator based quan- 
tum fluctuations on superconductivity. This began a new 
activity of gauge theory approach to strongly correlated 
electrons systems, and in particularly cuprate supercon- 
ductors. Some key ideas from lattice gauge theory was 
very effectively used to understand the strong correlation 
problems in Mott insulators and doped Mott insulators. 

This theory observed that the low energy Hamiltonian 
of Mott insulator has a local U(l) gauge symmetry, when 
expressed in terms of the underlying electron variables. 
This is manifest in the electron representation (equation 
5). A local (site dependent) gauge transformation 

4 - e'^'cL (22) 

leaves the spin Hamiltonian invariant. Because, blj trans- 
forms as 

bh ^ (23) 
and leaves superexchange hamiltonian (equation 5) 
—JJ2(ij) ^Ij^ij invariant. An important consequence of 
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this is a possibility of an emergent U(l) dynamical gauge 
field on the links. This link gauge field is a dynamically 
fluctuating complex field, Ay(t), connecting neighboring 
sites ij. Around the same time, while studying Anderson 
lattice and heavy fermion problems, Noga|64l] introduced 
such dynamically generated link variables. This insight 
was very useful for us and we used this effectively in our 
Mott insulator and doped Mott insulator problem. 

The dynamically generated field Ay was shown in 
path integral treatment of the problem, using a Hubbard- 
Stratanovic method, which reduces a quartic two body 
interaction term to a quadratic interaction term with dy- 
namically fluctuating pair field Ay(t). Path integral ex- 
pression for the partition function of the t-J model is: 

iar 

X Si d-r[-tT.(.,) 4(r)c^,,(r) + h.c.-JE(.,) b\^{r)b^ 

Here c|^(t)'s are Grossman variables; double occupancy 
constraint is formally expressed through a delta function 
along with the measure. Using a Gaussian identity, e" — 
e~'^^ ~'^^°-^dx we re-express the partition function 

as 

Z ^ j nrfA:,(T)dA,,(r)e-/o'''^-^5:,,>|A.,(r)|^ 

iar 

The auxiliary variable Ay (r) is the probability ampli- 
tude of finding a singlet bond between two electrons at 
sites ij. Now one can formally integrate over the fermi 
fields to get an effective action in terms of Ay (T)'s: 

J Yl DA*^{t)DA,j iT)efo rfrSe„[A.,(r)] (26) 

Without performing an explicit integration one can 
make some very general and useful statement about the 
symmetry of the effective action Scff[Aij (r)] at half filling 
and away from half filling. 

At half filling the Heisenberg model exhibits a local 
U(l) gauge invariance. This implies that 5'off [Aij(T)] is 
invariant under the local U(l) gauge transformation: 



A* 



e'^'A*^.e'"^- 



(27) 



fashion and obtained the following lattice action: 

'S'Mott ~ flO ^ ^ij^jkA-liAii + H.c. 

(ijkl) 



(28) 



This action exhibits a manifest local gauge invariance on 
the lattice. First term of the above action represents 
resonance of two valence bonds in an elementary square 
plaquette {ijkl). This is a large term, that is act ually 
proportional to J, as was noted later by Anderson 65 1. 
This term stabilizes resonating singlets in the ground 
state. When we focus on the soft phase variable 6',j, 
we can approximate the hard variable, |A.y| « Ao| and 
Ay « |Ao|e'^'J . In terms of the phase variable the U(l) 
action has the form: 



Smo 



clAol 



cos{9i. 

(ijkl) 



fjfe 



Jkl 



9h) (29) 



We evaluated the static part of the effective action, taking 
care of the double occupancy constraint in a mean field 



This paper summarized results of an analysis, using 
Villain approximation, and a plaquette integer vari- 
able (similar to vorticity in XY model in 2D) was dis- 
cussed. This integer plaquette variable was called 'mag- 
netic charges', using analogy to U(l) lattice gauge theory 
in 2+1 dimensions. These magnetic charges were iden- 
tified as key 'topological excitations' that will describe 
evolution of spin liquid state as a function of tempera- 
ture. The magnetic charges or vorticity, is a novel topo- 
logical excitation that arises in a disordered spin liquid 
state in 2D, the very first one as far as the authors knew. 
Kivelson-Rokhsar-Sethna's holon or spinon was shown to 
arise when Agj = 0, when they are located at the origin. 

It was further suggested that as the gauge field 
is dynamically generated, by interacting fermions, the 
magnetic charges might induce fermions, through a 
possible topological term. This will be a 'neutral 
fermion' (spinon) excitation. In the modern language 
it will be called a fermion-flux composite. Very soon 
Dzhyalozhinksi, Wiegman and Polvokov(55j suggested 
Hopf term as a topological term. It suggested a possibil- 
ity of statistics transmutation and a debate started about 
existence of Hopf term in 2D spin- ^ Heisernberg anti- 
ferromagnet. In parallel, Marston[6a| and also Zou 67 
suggested a Chern-Simon term in the U(l) gauge the- 
ory. Thus the gauge theory paper planted seed for the 
discussion of Chern-Simons field theory to describe 2D 
quantum spin systems. 

Soon the magnetic charge was called by a more ap- 
propriate name, 'magnetic flux' and Affleck and Marson 
found the 7r-flux RVB mean field solution in the BZA 
theory. This state respects parity and time reversal, as 
+7r and — tt can not be distinguished quantum mechani- 
cally (Bohm-Aharanov effect). A chiral spin liquid state 
that was suggested by Kalmayer and Laughlin|57|, in a 
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triangular lattice antiferromagnet was shown to be re- 
lated to an RVB mean field solution with condensed ^ 
flux, by Feng and LeejGSj. In more recent works the 
magnetic flux was christened as 'visons' in the works of 
Senthil and Fisher [s^l- Very recent work by GBj39], fol- 
lowing an early work of us with Anderson, John, Doucot 
and Liang [tO], finds that the well known skyrmion so- 
lution of 2D Heisenberg antiferromanget represents two 
unbound spinous that carry quarter magnetic flux each, 
and showing an important result that spinous are decon- 
fined semions, but above an energy gap. 

The magnetic charge or flux, that came in a natu- 
ral fashion in the RVB theory was described by Wen- 
Wilczek-Zeeji^ in terms of spin chirality variables ex- 
pressed, in terms of Pauli spin operators as: 



(30) 



Another important consequence of U(l) gauge field de- 
scription in a lattice was an application of Elitzur's 
theorem 7l| . Elitzur theorem states that in a pure lattice 



gauge theory (without matter) with local gauge invari- 
ance, the local gauge symmetry can not be spontaneously 
broken. Thus Elitzur theorem automatically precludes 
spontaneous symmetry breaking of the U(f ) gauge field 
in the Mott insulating state. It means absence of super- 
conductivity in a Mott insulator. This paper also pointed 
out possibilities of confined and deconfined phases of U(l) 
gauge theory, through the area and power law behav- 
ior of Wilson loop like loop correlation functions such as 
W{C) — (Tie ^ij^jk^ki---^pi) ■ In the deconfined phase the 
neutral fermions carrying spin-i moments were suggested 
to be present as low energy excitations. 

Next important question is the consequence of dop- 
ing and how it will bring about superconductivity in 
the gauge field description. This paper showed that the 
doped Mott insulator described by the t-J model does 
not have local U(l) symmetry of the parent Mott insu- 
lator. That is, in the presence of a non zero doping x 
0, the effective action ^off [Ay (t)] is invariant oly under 
a global U(l) gauge transformation, Ay — > e*^^Ay. 

We evaluated the effective action, using the same ap- 
proximation as before and obtained, 

S ^ ao^A*^.A,fcA*;A„; + ff.c. 

i'jkl 

+ xai'^A*jAjk + H.c. 

{ijk 

+ 6^|Ap + c^|Ay|4 
{ij) (*i> 



(31) 



This had a remarkable consequence. As the second term 
proportional to x, the dopant density, removes the lo- 
cal gauge invariance, total action has only a global U(l) 
gauge symmetry. The first term representing a plaque- 
tte resonance of the valence bonds is the memory of the 



Mott insulator in this approach. An immediate conse- 
quence was that Elitzur's theorem is no more applicable 
now. In principle superconductivity is possible. Action 
(equation 31) is the simplest lattice Ginzburg Landau 
action for RVB superconductors. The major aspect of 
RVB appears from the plaquette resonance term, which 
fights against long range order. During 1988, Nakamura 
and Matsui[22| used the above lattice action and did a 
complete numerical evaluation of the partition function 
(going beyond saddle point approximation) and found 
reasonable superconducting Tc close to doping. This was 
an important numerical proof that 2D superconductivity 
survives gauge field fiuctuations at optimal dopings. 

Soon after the gauge theory paper, MuUer-Hartman's 
group [ll] calculated the coefficient of the static lattice 
RVB-GL action, taking care of double occupancy con- 
straint more accurately and obtained a detailed phase 
diagram in the x-T plane. 

The effective action found the gauge theory paper, like 
the BZA method, has the right physics as far as sym- 
metry of the superconducting order parameter is con- 
cerned. The coefficient of the 'gradient term' in the lat- 
tice xai J2{ijk ^ij^jk+H.c. had the right sign ( xai > 0). 
And minimization of the above auction automatically 
leads to a d^^-y^ solution, very similar to Kotliar's d- 
wave solution. 

The RVB-GL action above does not have fermions ex- 
plicitly, as they have been integrated out completely. 
This is not correct for a d^^-y^ superconductor as it has 
nodal fermion quasi particles. So the low energy Bogoli- 
ubov quasi particles and their coupling to A^ should be 
present as part of the action. This is easy to incorpo- 
rate either phenomenologically or microscopically. Many 
authors such as, Affleck, Marston, Matsui, Nakamura, 
Wiegman, loffe, Larkin, Patrick Lee, Wen, Wilczek, Zee, 
Nagaosa, Read, Sachdev, Balents, Nayak, Fisher, Senthil, 
Dung-Hai Lee, Tesanovic, Franz and others have con- 
tributed to the elaboration of these fundamental ideas, 
most of them using continuum action. 

However, it is important to point out that it is difficult 
to incorporate the plaquette term, which keeps the mem- 
ory of the Mott insulator, in a continuum approximation. 
To this extent, the RVB-GL theory on a lattice remains 
unexplored. We find that because of the plaquette term, 
which distinguishes it from a standard superconductor, 
many interesting consequences could occur; for example, 
Andreev bound states at the vortex core and nature of 
impurity states induced by Zn and Ni substitution at 
copper site. In a recent work Muthukumar and Weng[73l| 
have developed an RVB-GL theory, starting from a slave 
fermion approach, and studied the physics of spinon vor- 
tices and properties of electromagnetic vortex core. 

Completing the BZA Programme 

1. d-wave Solution in BZA theory 
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As mentioned earlier, the available experimental data 
at the beginning of 1987 was not inconsistent with ab- 
sence of long range magnetic order in the Mott insulating 
end. Further, in the doped Mott insulator, experiments 
continued to show a small amount of electronic linear 
specific heat at low temperatures in the superconducting 
region. Both these results gave PWA a confidence that 
the quantum spin liquid with a pseudo fermi surface at 
the Mott insulator continues to become a superconductor 
with extended-s symmetry, and retains neutral fermionic 
excitations with a pseudo fermi surface. Overwhelmed 
by this confidence we were sailing happily in a psuedo 
fermi sea. 

During this period, intrigued by the structure of BZA 
mean field theory and the gauge theory^ approach, Af- 
fieck and Marston 2]| found a tt RVB fiux mean field solu- 



tion for the Mott insulator and Kotliar[17| the dx2_y2 so- 
lution for the superconductor, within the BZA approach. 
Affleck and Marston presented their result as an exact re- 
sult in a large N limit of a generalized Heisenberg model 
in 2D. 

As mentioned earlier, mean field energy of d-wave be- 
ing lower than extended s-wave is seen as follows. In BZA 
theory, the pair potential, for pair scattering, is given by 
-J7(k - k') = --J[cosik^ - k'J + cos{ky - k'y)]. This 
potential, which is attractive for small momentum trans- 
fer k — k' , changes sign and becomes repulsive for large 
momentum transfer k — k' (±7r, ±7r). The Ax^-yi so- 
lution, which changes sign as we move along the fermi 
surface in k-space takes advantage of this and becomes a 
lower energy state. 

Even after Kotliar's solution, there was a reluctance 
from Princeton group to accept d-wave solution. In ad- 
dition to the then existing phenomenological support for 
extended-s solution, we were not sure if the BZA d- 
wave solution will continue to have lower energy than 
extended-s solution, after Gutzwiller projection. Later 
Gutzwiller projected d-wave BCS solution was shown to 
have lower energy using numerical methods. 

The beauty of Affleck-Marson and Kotliar's solution 
was that both had nodal quasi particles. Various exper- 
imental results that followed later demostrated clearly 
d-wave superconductivity for hole doped cuprates. Need- 
less to say that an unconventional order parameter 
such as d-wave has its own profound consequences for 
cuprates, as has been seen both in the experimental and 
theoretical fronts in the last two decades. It is fair to say 
that RVB theory contains in its bag all these fascinating 
possibilities and perhaps even more. A recent excitement 
is with respect to superconductivity in Naa;Co02-yH2 0. 
We has developed an RVB theory for this system[7J|, 
where d -|- id, another unconventional order parameter 
that violates parity and time reversal symmetry is pre- 
dicted. 

2. Improving BZA mean field Hamiltonian 



As discussed earlier, the hard problem of Gutzwiller 
projection was replaced by an ansatz in the variational 
analysis. We replaced the hopping parameter t by tx. 
This tells us that an electron can hop to a neighboring site 
only when it is empty. The probability of it being empty 
is X, the doping density. This is equivalent to replacing 
bare electron mass by a renormalized mass, which was 
already familiar to us from the work of Brinkman and 
Rice in the context of Mott insulator to metal transition 
in Hubbard model at half filling. 

While evaluating expectation values involving 
Gutzwiller projection, an approximate method that 
takes care of some incoherent aspects of projection, was 
developed by Gutzwiller, in his earlier study of Hubbard 
model. It involved certain combinatorics. This method 
was already successfully used by Rice,Joynt, Shiba, 
Ogatta, Volhardt, Anderson and others in dealing with 
issues of heavy fermions and also Hubbard modeling of 
solid '^He. This method was adapted by Zhang, Gros, 
Rice and Shiba^2Q] to improve the BZA renormalization 
process. Their improvement is schematically shown 
below. 

PGiHt+Hs)PG ^-^^ (xHt+Hs) (gtHt+gsHs) 

(32) 

with doping dependent renormalization parameters gt = 
and gx — jYf^p^- But for the two new renormal- 
ized parameters the rest of the self consistent theory was 
identical to BZA theory. Further, physical correlation 
functions such as anomalous amplitudes (blj), got an ap- 
propriate renormalization factor gt times mean field am- 
plitude : 



(4) 



.9t(&I,)o 



(33) 



This is in spirit similar to BZA theory, where it is stated 
that fraction of singlet bonds that are charged are x and 
use this to calculate superconducting Tc . It should be 
also pointed out that the renormalization parameters gt 
and gs are also not variational parameters, very much 
like in BZA theory. Here also one gets a gap equation for 
the order parameter A and p, which can be solved self 
consistently. The merit of the improvement suggested by 
ZGRS is that the results obtained by a mean field anaylsis 
is quantitatively more accurate. This is elaborated in the 
Plain Vanila paper as well as more recent works of Zhang 
et al and Muthukumar-Gros et al. 

Generalization of the above BZA program to fi- 
nite temperature has problems. This problem is re- 
lated to variational states that are orthogonal before 
Gutzwiller projection becoming non-orthogonal after 
Gutzwiller projection and decrease in density of rele- 
vant Hilbert space introduced by double occupancy con- 
straint. Anderson |ilS| has recently offered an analysis, 
where he introduces the notion of 'spin-charge locking' 
and introduces a generalized BCS type of formalism, 
where separate Anderson-Nambu spinors are introduced 
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to take care of neutral spin-pairing and electron pairing. 
At very high temperatures they are decoupled and in the 
superconducting state they are locked. Spin-gap phase 
represents a progressive locking. 

3. Holon Condensation &c Slave Boson theory - 
a support for BZA program 



Soon after Anderson's paper and BZA theory, an in- 
sightful paper by Kivelson, Rokhsar and Sethna[l3| of- 
fered a theory where a doped charge enters as holon, 
a spin zero soliton with charge -l-e. They act like 
bosons and bose condense leading to superconductivity. 
It seemed to express, in a formally correct fashion, An- 
derson and BZA's result of viewing the charged valence 
bond as a boson and their bose condensation. Zou and 
AndersonQ adapted the slave particle method, devel- 
oped earlier by Barnes, Coleman, Read, Newns, Kotliar 
and Ruckenstein, to the t-J model. This began a series 
of investigations starting with the work of Kotliar-Liu, 
gauge theory by Wiegman[23j| and the present author 24 1, 
and detailed mean field analysis by Fukuyama[l^ and 
collaborators and several other investigations. 

The advantage of slave particle method over the vari- 
ational approach is that one can introduce dynamically 
generated gauge fields and go beyond mean field theo- 
ries, without constrained by variational RVB states and 
also address finite temperature problems. However they 
are technically hard, as it is clear in the works of many 
authors that followed. BZA mean field theory seems to 
capture the essence of superconductivity phenomena in 
the t-J model. 

In some limit the slave boson analysis gives the same 
result as the improved BZA Hamiltonian by Zhang et 
al. During this development, we also showed[6| that the 
holon condensation is not actually a charge 2e rather 
charge e condensation in view of the double occupancy 
constraint. In other words, electron pairs get effectively 
delocalized into a zero momentum condensed state. Or 
holon is book keeping device for a correlated fluctuation 
and delocalization of charged valence bonds, as far as 
superconductivity is concerned. 



4. Lurking dangers outside Optimal Doping 

As RVB theory was being developed, there was a con- 
scious effort to focus on region around optimal doping. 
RVB superconductivity is at its best and one can hope 
to understand it better here, than elsewhere. In fact, t- 
J model is a reasonable model for real cuprates, 
when we have either an isolated hole in a Mott 
insulator or a finite density of holes at optimal 
doping ! Any thing in between is complicated because 
of unscreened long range coulomb interactions among 
charge careers, strong coupling to phonons and disor- 
der effects from off plane dopants. It could easily lead 



to nano scale phase separation into Mott insulating and 
optimally doped regions, self trapping etc. 

In reality it turned out to be even more complicated. A 
doped Mott insulator also supports other quasi or psuedo 
long range orders such as fluctuating charge and spin or- 
ders and chiral orders, involving circulating spin currents 
and charge currents. Interactions not contained in the 
t-J model, seem to encourage these competing orders, 
outside the superconducting dome. They are fascinat- 
ing quantum and classical condensed matter problems. 
However, they are not crucial in our foremost goal of un- 
derstanding high Tc superconductivity deeply. There are 
good phenomenological and theoretical evidences that 
they are not the root cause of high Tc superconductivity. 
In fact, they are competitors [sij 

The same is true of the so called pseudo gap phase, out- 
side the dome. Experimental evidences point to presence 
of enhanced correlations corresponding to the cometing 
order. 

There is a clear indication that the physics and inter- 
actions contained in the t-J model 2D is alone able to 
support high Tc superconductivity, without need for an 
external help or a catalytic agent. 

Mean field theory works remarkably well 
at optimal doping 

RVB mean field theory and gauge theory approach, 
predicted a robust superconducting state at optimal dop- 
ing with a large Tc and supported, Anderson's proposal. 
These papers also pointed out that in the superconduct- 
ing region, mean delocalization energy per particle (~ xt) 
and energy of superexchange (J) are comparable. The su- 
perconducting condensation energy is a finite fraction of 
J or xt. Now we know that ^ ISOOii' for cuprates and 
superconducting condensation energy is large and so are 
superconducting Tc 's. This puts cuprates on a different 
region in the Uemura plot. In fact, Uemura plot itself 
was inspired by RVB type of idea of bose condensation 
of charged valence bond in cuprates. 

Fortunately for the above theories, new cuprate family 
members were discovered, where superconducting Tc 's 
soared to new heights. The Hg and Tl based single layer 
superconductors reached a Tc in the range of 95 K. At 
optimal doping, Tc is large compared to its one layer 
counterpart LSCO or one layer BISCO material. I have 
argued elsewhere, that a single layer material has a large 
intrinsic superconducting Tc w 120K. Competing or- 
ders (charge and spin order tendencies, often helped by 
octahedral rotation or distortions) steal away the super- 
conducting condensation energy, making superconduct- 
ing Tc among one layer materials to swing from 5 K to 
95 K at optimal doping. 

The highest Tc is in the Tl based cuprates, where Tc is 
as large as 163 K, under large external pressures. These 
are remarkable experimental support to the BZA and 
gauge theory approximations. 
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Several theoretical attempts readily show enhanced 
singlet correlations in the ground state, a prerequisite for 
long-range singlet superconductivity. Variational Monte 
carlo analysis of RVB superconducting wave functions 
have given useful results, where quantitative comparisons 
have been made with some experimental quantities. Over 
years quantum Motne carlo methods have given encour- 
aging results. Other approaches, such as real space as 
well as k-space cluster DMFT have given very meaning- 
ful results, consistent with BZA RVB mean field theory. 

One of the best support for BZA program was provided 
by a serai phenomenological theory due to Emery and 
Kivelson|62l|. which focused on superconducting transi- 
tion temperature. Kivelson and Emery assumed, con- 
sistent with RVB theory and the existing phenomenol- 
ogy that superconducting phenomena at low doping is 
dominated by phase fluctuations rather than amplitude 
fluctuations. This leads to a simple expression for super- 
conducting Tc in terms of measurable quantities. In par- 
ticular energy associated with spatial variation of phase 
is expressed as: 

H « / {mr)fdr (34) 

Here m* is the effective mass of cooper pair and n^. (0) is 
the superfluid density at zero temperatures. Then one 
can use the standard KT formula to get an expression 
for superconducting Tc as 

^^^^ " 

This expression is remarkably similar to the expression 
given in BZA paper and the later paper with Zou and Hsu 
(equation 20 and 21). Emery and Kivelson went further 
and made detailed comparison with existing phenomenol- 
ogy, including how Tc increases with inter layer coupling 
etc. The remarkable result is that, as in the BZA result, 
here also large U or J does not explicitly appear in the fi- 
nal result for ■ Superexchange, a consequence of large 
U, simply provides a large energy scale and an umbrella 
below which fairly stable singlets can delocalize and pro- 
duce a superconducting state, along the lines of BEG. As 
thermally produced nodal quasi particles will interfere 
with superconductivity, in a way different from the fully 
gaped superconductors, certain corrections were neces- 
sary to the above analysis, as indicated in reference [33j|. 

Even though no rigorous theorem exists proving exis- 
tence of a finite Tc superconductivity in large U repul- 
sive Hubbard model or t-J model in 2D, it is increasingly 
becoming clear that such a theorem is likely to exist. 
However, one should remember that spin-i systems are 
notoriously hard. For example, a rigorous demonstration 
of long range antiferromagnetic order, in the undoped 2D 
square lattice Heisenberg antiferromagnet does not exist 
even now. Will doping make it any simpler ? 



As mentioned earlier, in the heat of developments dur- 
ing the beginning of 1987, a phase diagram was conjec- 
tured in the paper by Anderson, Hsu, Zou and the present 
author, based on BZA theory (figure 3). This phase di- 
agram was conjectured well before experimental phase 
diagram emerged. It is interesting that the experimental 
phase diagram had an excellent overal form as conjec- 
tured by the theory. This again indicates that the fluc- 
tuation effects do not change the qualitative prediction 
of the mean field theory. 

BZA paper also pointed out that ' ... in the low dop- 
ing concentration limit (x << 1) phase fluctuations play 
an important role and the mean field theory fails. At 
this limit, Tc and gap are governed by the large phase 
fluctuations. ... But for large enough x (> 5%) we ex- 
pect that the mean field theory works.' The main reason 
for the conjecture was that, in order to be able to con- 
struct a phase coherent superconducting state, two most 
important requirements are i) well developed short range 
singlet correlations and ii) finite long wavelength charge 
compressibility. Singlet correlation is there in plenty in 
the doped Mott insulator, because super exchange con- 
tinues to be present for a range of doping. A doping of 
5 % gives sufficient delocalization energy to escape local- 
ization effects due to disorder or hole-hole electrostatic 
interaction and attain a finite charge compressibility. 

Spin-1 collective mode at (tt, tt) first seen in neutron 
scattering [t^, is one of the key unconventional feature of 
cuprate superconductivity. As far as I know no such spin- 
1 collective mode is seen in any conventional supercon- 
ductor. Moreover, the frequency of this mode scales lin- 
early with superconducting Tc ■ In fact, when Tc becomes 
zero this mode becomes the Goldstone mode of the an- 
tiferromagnetic order, that is supported by a Mott in- 
sulator through superexchange. Thus spin-1 collective 
mode is a memory of the Mott insulator. It is a 
simple manifestation of the tightly bound bond singlets 
that makeup the superconducting state: these singlets 
are softer at (tt, tt) in k-space. This reveals a dynamic 
antiferromagnetic correlation at (7r,7r), at energy scales 
large compared to the superconducting energy scale. The 
simplest theory that explains the spin-1 collective mode 
at (tt, tt), in a natural fashion at (tt, tt) is an RPA collec- 
tive mode analysis built on a BZA type of d-wave mean 
field background, within a t-J model. This is another 
support for the mean field theory. 

In a recent paper Anderson and collaborators |48|] have 
made comparison of the physical properties calculated 
within ZGRS approach, variational Monte-Garlo analy- 
sis and experimental results. The agreement between the 
three are very good, underscoring the validity of the orig- 
inal BZA approach and the assumptions therein. 

Another reason why a simple BGS type of theory works 
well at optimal doping is the following. For zero and 
very small doping the strongly correlated electronic sys- 
tem has strongly localized electrons and there is no re- 
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semblance of the ground state to a metallic fermi sea. 
For very high doping we clearly have a degenerate fermi 
gas; at these heavy dopings superexchange does not take 
place, as electrons have a large mean velocity and do not 
have time for superexchange interaction with a neighbor 
in real space. At optimal dopings a reasonable fermi sea 
is formed at the same time superexchange also survives. 
It is this combination which makes the situation very 
similar to a BCS theory with superexchange being the 
glue. This is the reason BZA renormalized Hamiltonian 
and the improved version by ZGRS work well. 

Is Spin Fluctuation a Glue ? 

The BZA program that was initiated in 1987 and got 
completed shortly, is capable of answering many ques- 
tions in the superconducting region. Surprisingly such 
work began only recently. Partly because, we were either 
ignorant or ignored key old developments. Too many 
questions, not necessarily related to the major debate, 
arose from a wealth of experimental results. Being a com- 
plex system, cause and effects get mixed up sometimes. 
We will use the spin fluctuation scenario to illustrate our 
point. 

Historically, Hirsch's numerical analysis of the single 
bond Hubbard model, brought out the beautiful possi- 
bility of an extended-s wave superconducting correlation. 
This and other developemens started focusing on this 
phenomena in terms of diagrammatics and exchange of 
spin fluctuation bubbles. 

Let us look at experimental facts. At optimal doping, a 
spin-1 resonance emerges at (tt, tt), as a sharp mode, but 
only in the superconducting state. In the normal state 
it becomes very broad and disappears. This is a col- 
lective mode, characteristic of the RVB superconducting 
state. Further, the energy of this mode scales linearly 
with superconducting Tc . Formally, as the supercon- 
ducting Tc vanishes, this mode becomes a spin wave or 
Goldstone mode. Further, in addition to the spin-1 res- 
onance there is a broad background of spin fluctuation 
activities, as seen in S(q,a;). 

RVB theory takes the point of view that spin flucua- 
tion contribution to S(q,CL') around (tTjTt) is a result of a 
built up of singlet correlation through superexchang pro- 
cesses. These fluctuations that occur in a broad range in 
momentum and energy space are far from any coherent 
modes. Instead of talking about these dissipative spec- 
trum of spin fluctuation activities, RVB theory directly 
focuses on superexchange J and treats it as a glue. The 
only coherent mode is the spin-1 resonance. As it has 
been pointed out, it has a very small spectral weight and 
further it occurs only below Tc . So it can not be a glue 
either. 

All these phenomena are manifest in diagrammatic 
spin fluctuation calculations. Electron has a strong fre- 
quency and momentum dependent and large normal and 



anomalous self energies. They depend on each other self 
consistently. The large normal part of the self energy 
makes the quasi particles spectral function very broad, 
consistent with ARPES experiments. But this approach 
misses the built up of singlet correlations, as a result 
there is no way we can even approach under doped re- 
gion, leave alone the Mott insulator region by this ap- 
proach. There is no natural and simple way of getting 
supcrfluid density being proportional to doping at small 

X. 

There is a missing logic. While one may get some sat- 
isfactory answers, a complete picture, a claimed strength 
of RVB theory, is absent. 

Conclusion 

In condensed matter physics we attempt to synthesize 
new laws, notions and introduce new reference states or 
phases in demystifying properties of complex materials. 
There are several idealized reference states for describing 
a variety of quantum phenomena in solids: free electron 
gas, ideal bose gas, harmonic phonons, fermi liquids, Lut- 
tinger liquids, BCS paired fermi sea, etc. The wealth of 
quantum condensed matter phenomena force us to intro- 
duce new reference states occasionally. 

PWA's RVB proposal in 1987 and subsequent devel- 
opments in the last 20 years illustrates a struggle to in- 
troduce Mott insulator as a reference phase to describe 
the unusual high Tc superconductivity and a variety of 
related anomalies Q in the family of cuprates. From su- 
perconductivity point of view, it is a serious attempt to 
find an alternative to phonon mediated superconductiv- 
ity, compelled by experiments. 

Looking back, it has been a worthwhile struggle, and 
RVB theory has silently entered the subconscious of the 
condensed matter mind. La2Cu04 has become a text 
book Mott insulator. A one band model, with strong 
correlation, namely t-J or large U Hubbard model is ac- 
cepted as a minimal model to describe low energy physics 
of cuprates. Deep consequences of projection in these 
models are also getting accepted, but slowly. The U(l) 
RVB gauge field theory Q has been nurtured and devel- 
oped by several authors in commendable ways. Old ideas 
from RVB are being rediscovered. The pseudo fermi sur- 
face and neutral fermions sugg ested PWAd and shown 
to be a possibility in a microscopic theory in BZA, is 
bein g r ealized in certain organic Mott insulator, ET- 
salts[76|. RVB mechanism is being successfully applied to 
understand superconductivity in organics 1781 1. A new 
superconductor Na:ECo02- yH20, is likely a long awaited 
doped spin-i Mott insulator on a triangular lattice [t^. 
There is a strong indication for RVB superconductivity in 
boron doped diamond (7!^]. through superexchange effects 
in an impurity band Mott insulator. 

In the field of superconductivity, before cuprates ap- 
peared in the scene, one was used to the luxury of 
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a beautiful and powerful BCS theory, that is so use- 
ful in understanding majority of elemental supercon- 
ductors. Formalisms such as the Eliashberg theory, a 
microscopic approach and other related developments 
have been extremely helpful in understanding many 
experimentally measured properties, including tunnel- 
ing spectra, a^Fiu;) etc., accurately. Ginzburg-Landau 
phenomenological theory gets a microscopic meaning 
through Nambu-Gorkov formalism. Space and time de- 
pendent superconductivity phenomena and quasi particle 
dynamics is understood through Bogoliubov-de Gennes 
theories. There are new phenomena that came as predic- 
tion after BCS theory: Josephson effects, Andreev reflec- 
tion etc. All is well with old superconductivity. Once we 
agree to live with a few parameters of less microscopic 
origin, even some bad actors such as A15 and Chevral 
phase compounds seem to yield to BCS theory. 

The situation with cuprates and many new materi- 
als, suspected to be RVB superconductors are different. 
We mentioned about the lurking dangers out side op- 
timal doping from competing phases and extra interac- 
tions that could encourage their growth. This makes a 
theory based on t-J model of some what limited validity 
! Much care should be taken to get all the low energy 
physics from the t-J model. For example, determination 
of superconducting Tc for a given system, say LSCO is 
wrought with complications in the way we explain below. 

In BCS theory we used to worry about isotope shift of 
a fraction of a degree in Tc . In cuprates, at optimal dop- 
ing, in one layer materials, the Tc varies between 5 K to 
95 K between Bismuth and Thallium one layer materials. 
All these materials have very similar normal state prop- 
erties, such as the coefficient of linear resistivity, modulo 
some material purity complications. Such a large varia- 
tion indicates that superconductivity is not alone. There 
are other factors and competitors that are at work. Sin- 
gle layer superconducting Tc may even get enhanced by 
an interlayer pair tunnelling phenomena[8| in bilayer sys- 
tems, an entirely new additional contribution that owes 
its origin to anomalous normal state. 

So the world of cuprate superconductors are different 
and complex. First one has to answer some generic ques- 
tions as accurately as possible for the 2 dimensional t-J 
model. For example, what is the x dependence Tc near 
optimal doping. What is the maximum Tc at optimal 
doping for a range of t and t' and J that is relevant for 
cuprates. How Tc gets modified in a bilayer or multilayer 
system for a given tj_. 

It is not meaningful to pick up LSCO and try to un- 
derstand superconducting Tc from t-J model. We will be 
tempted to find the correct band parameters t, t', J and 
do an accurate RVB mean field theory. It is very clear 
that what controls the experimentally seen reduced Tc is 
some phenomena outside t-J modeling or band structure 
effects. In fact, in one experiment [80] an epitaxial strain 
increases the superconducting Tc of LSCO thin film from 



25 to 50 K, without any manifest increase in doping. Sim- 
ple estimate shows that epitaxial strain can not give suf- 
ficient change in t's to cause such a large effect. The 
fundamental reason seems to be that epitaxial strain is 
encouraging certain atomic scale lattice distortion that 
encourages charge order formation self consistently. 

There is some new physics 81| , namely a competing 
phase, which is stealing away superconducting condensa- 
tion energy. As BISCO and Tl one layer materials are 
affected in different ways, actual theoretical prediction 
of ground state gap parameter, will require additional in- 
puts, either phenomenological or microscopic. One of the 
important question is how easily the Cu-Oxygen octahe- 
dron, or pyramid or square planar complex respond and 
get distorted or rotate in their different environments, to 
either charge or spin or valence bond localizations and 
ordering tendencies. Then there will be feed back and 
growth of competing orders, at the expense of supercon- 
ductivity. 

One of the strategies will be to pick the best super- 
conductor in the single layer family for a deeper un- 
derstanding of mechanism of superconductivity. From 
this point of view our first superconductor namely doped 
La2Cu04 is a bad system to study ! Superconducting 
Tc never goes beyond 30's at optimal doping. Do we 
have to abandon two decades of experimental efforts ? 
Perhaps not. Then an important question is why doped 
La2Cu04 , or other low Tc members have never attained 
their full potential (maximum Tc ) that Tl one layer sys- 
tem has reached. These are new questions that has no 
parallel in conventional superconductors. 

At the end it may be possible to accommodate effects 
of competing orders phenomenologically by modifying t', 
J etc. But that might miss important physics. It is 
important to recognize and treat competing dynamical 
processes in their own right. This of course makes the 
problem hard. 

So we realize after 20 years of cuprate study that we 
are in a different situation. The nature of questions asked 
should be different and method of analysis will be differ- 
ent. Priorities will be different. For example one of the 
important question will be how to reach the full potential 
or maximum superconducting Tc in the one layer family ? 
As we said earlier, the one layer superconducting Tc can 
be as large as 95 K. However, in multi layers we have a 
superconducting Tc of 163 K. Is it a consequence of coop- 
eration from interlayer pair tunneling phenomena within 
the multi layers, or some structural rigidity that discour- 
ages competing orders, there by making a single layer 
realize new heights ? 

We have focused on only one thermodynamic property. 
We have similar question about the energy of the spin-1 
collective mode and a variety of other key physical prop- 
erties. 

There are several important experimental issues in the 
superconducting state: temperature dependence of order 
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parameter, nodal quasi particle dynamics, as revealed by 
magnetic resonance studies, S(q,a;) from neutron scat- 
tering, spectral functions of quasi particles in ARPES, 
detailed space and energy dependent STM study of local 
electronic density of states and gap functions, quasi par- 
ticle interference effects, structure of vortex core, states 
inside the vortex core, bound states around impurities 
such as Zn, Ni and their different and unexpected ef- 
fects as a function of temperature, thermal conductivity 
anomalies etc. etc. 

Suddenly one finds oneself in the midst of a flood of 
questions and real complications, so different from ele- 
mental superconductors. As we mentioned earlier it is a 
meeting place of quantum magnetism and superconduc- 
tivity, Mott insulator and fermi sea. Many things other 
than superconductivity are taking place. It is unfair to 
say we do not understand high Tc superconductivity. We 
understand it too well, so we fear it and tread carefully. 

Anomalous normal state phenomena takes us to a dif- 
ferent world all together. Experimentally, it is a clear 
case for a non fermi liquid in 2D. The many body theory 
we have developed, to address these challenges, are at the 
beginning stage. PWA has attacked this problem from 
different angles : i) a non- vanishing phase shift at 2ki?, in 
the forward scattering spin singlet channel, in 2D Hub- 
bard model ii) tomographic Luttinger liquid, asymptotic 
Bethe ansatz in 2D iii) Anderson-Khveschekno's anoma- 
lous commutation relations, iv) spin charge decoupling 
and two relaxation times on the fermi surface v) asymme- 
try in single particle tunneling arising from the key pro- 
jection in a t-J model and vi) very recently an orthoganal- 
ity catastrophe inherent in the projected t-J model etc. 

Similarly, in defect free underdoped cuprates super- 
conductivity seems to vanish because of quantum fluc- 
tuations or unscreened coulomb interactions, leaving a 
metallic ground state. This possibility that was sug- 
gested earlier with a pseudo fermi surface and later 
with nodal quasi particles [s^]. The one with nodal exci- 
tations seems to be gaining experimental support [s^l in 
very pure YBCO. It will be a pristine RVB state, a ref- 
erence non fermi liquid state, that has not yielded to in- 
stabilities such as charge order or spin order or supercon- 
ductivity. P and T violating metallic ground states are 
likely to exist in organics and in cobalt oxide systems [3], 
near the Mott insulator end. 

So it is humbling to see the kind of problems and also 
richness that Bednorz and Miillers discovery and Ander- 
son's RVB proposal has created in the world of supercon- 
ductivity and strongly correlated electronic systems. At 
the same time it is heartening to see that BZA and related 
papers have answered the primary question of existence 
of high temperature superconductivity in cuprates in a 
microscopic theory, rather satisfactorily, and is ready to 
answer several old and new questions. 
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